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ABSTRACT 

We present ultra-deep Spitzer 70/ini observations of COODS-North (Creat Observatories Origins Deep 
Survey). For the first time, the turn-over in the 70/im Euchdcan-normahzed differential source counts 
is observed. We derive source counts down to a flux density of 1.2 mJy. From the measured source 
counts and fluctuation analysis, we estimate a power-law approximation of the faint 70/im source counts 
of dN/dS oc 5^^'^, consistent with that observed for the faint 2Afim sources. An extrapolation of 
the 70/im source counts to zero flux density implies a total extragalactic background light (EBL) of 
7.4 ± 1.9 nW m~^sr~^. The source counts above 1.2 mJy account for about 60% of the estimated EBL. 
From fluctuation analysis, we derive a photometric confusion level of CTc = 0.30 ± 0.15 mJy {q = 5) for 
the Spitzer 70/im band. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: high-redshift — infrared: 
galaxies 



1. INTRODUCTION 

Deep 24/im observations (Chary et al. 2004; Papovich 
et al. 2004; Fadda et al. 2006) have demonstrated the 
ability of the Multiband Imaging Photometer for Spitzer 
(MIPS, Rieke et al. 2004) to study the mid-infrared (mid- 
IR) properties of high-redshift galaxies (Yan et al. 2004; 
Le Floc'h et al. 2004, 2005; Perez- Gonzalez et al. 2005; 
Daddi et al. 2005; Papovich et al. 2006; Caputi et al. 
2006). The interpretation of the 24/im data are compli- 
[ cated by the presence of strong emission and absorption 
features (e.g., Armus et al. 2004) redshifted into the 24/im 
band. Observations at longer wavelengths, such as 70/im 
which is closer to the far-infrared peak of the spectral en- 
ergy distribution (SED) and is away from the strong mid- 
IR features, are crucial for constraining the infrared lumi- 
nosities and star-formation rates. 

The previous deep Guaranteed Time Observer (GTO) 
surveys did not achieve sufficient sensitivity at 70/xm to de- 
tect distant luminous infrared galaxies (LIRGs; 10^^ L© <^ 
Lir < 10^^ Lq), without stacking 70/im data for a large 
number of 24/^m-selected sources (Dole et al. 2006). Much 
deeper observations are needed at 70/im to individually de- 
tect the z ~ 1 LIRGs that account for the majority of the 
extragalactic background light (Elbaz et al. 2002; Lagache 
et al. 2004). In this letter, we present initial results for 
the deepest 70/im survey taken to date with Spitzer. 

2. OBSERVATIONS 



The ultra-deep 70/im observations of the northern field 
of the Great Observatories Origins Deep Survey (GOODS- 
N) were carried out in Cycle-1 of the General Observer 
(GO) program {Spitzer program 3325). The survey covers 
the central 10' x 10' of GOODS-N to a depth of 10.6 ks. 
The data were taken using small-field photometry mode 
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with lOMIPS-second data collection events (DCEs). The 
field was observed with an 8 position cluster map for each 
astronomical observational request (AOR). The observa- 
tions were repeated with 12 AORs taking 34.5 hours of 
observatory time in total. The mapping order and dither 
positions of the cluster positions within the AORs were 
varied to provide uniform coverage and data quality across 
the field. The data were embargoed until after the GTO 
proprietary period and were released to our team in 2005 
August. In addition to the GO data, we used the MIPS 
GTO data of GOODS-N {Spitzer program 81, Dole et al. 
2004a) . The GTO data were taken in slow scan mode with 
one degree scan legs and have an integration time of 600 s 
at TOfim, covering an area of 0.6 deg^. 

3. DATA REDUCTION 

The raw data were downloaded from the Spitzer Science 
Center (SSC) archive and were processed from scratch us- 
ing the offline Germanium Reprocessing Tools (GeRT, S13 
version 1.0). The instrumental artifacts in the basic cal- 
ibration data (BCDs) were removed adopting the filter- 
ing techniques used for the reduction of the extragalac- 
tic First Look Survey (xFLS, Prayer et al. 2006). The 
BCD pipeline processing and filtering procedures were op- 
timized for these deep photometry data. We adopted the 
updated S13 calibration, which assumes an absolute fiux 
calibration factor based on stellar SEDs of 702MJysr~^ 
per MIPS-70 data unit. We then multiplied the data 
by the color correction factor of 1.09 to place the data 
on a constant vfi, scale, which is also the appropriate 
color correction (within 2%) for a wide range of possi- 
ble galaxy SEDs (see SSC web pages for calibration and 
color-correction details). In comparison, the calibration 
correction adopted here is 3.4% larger than the calibra- 
tion adopted for the xFLS analysis (Prayer et al. 2006). 
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The filtering of tlie data is a crucial aspect in the pro- 
cessing. For the 70/im photometry mode, calibration stim- 
ulator (stim) flashes are used every 6 DCEs, and latents 
due to these stim flashes accumulate over time. To re- 
move stim flash latents and additional artifacts, we used a 
median column filter followed by a median high-pass time 
filter per pixel (with a filter width of 16 DCEs). The po- 
sitions of bright sources in the BCDs were fiagged so that 
the filtering corrections were not biased by the presence of 
sources. The median filtering techniques yield small off- 
sets from zero in the average level of the filtered-BCDs 
(fBCDs). These offsets correlate with the DCE position 
within the stim cycle and were removed by subtracting 
the median level from each fBCD. 

The data were coadded onto a sky grid with 4" pix- 
els using the SSC mosaicing and source extraction soft- 
ware (MOPEX, version 112505). Sources were extracted 
using MOPEX point source response function (PRF) fit- 
ting. In crowded regions, the 24yu.m positions (R. Chary 
et al., in preparation) were used for deblending. For op- 
timal source extraction, an accurate uncertainty image is 
needed. The uncertainty image was constructed by com- 
bining the noise per pixel based on repeated observations 
with the local spatial pixel-to-pixel dispersion after the ex- 
traction of bright sources. The average level of the uncer- 
tainty image was then scaled to match the average empir- 
ical point-source noise derived by making multiple aper- 
ture measurements at random locations throughout the 
residual mosaic after source extraction. The scale factor 
between the aperture derived point source noise (includ- 
ing the aperture correction) and the pixel surface bright- 
ness noise is iT(point source)/cr(4" pixels) = 10.9 ± 1.1 
[mJy/(MJy sr^"'^)]. The average point source noise for the 
ultra-deep area (after the extraction of sources) is 0.53 mJy 
(Icj). The applicability of the uncertainty image for point 
source extraction was verified by obtaining PRF fits with 
~ 1 for sources throughout the image. 

4. RESULTS 
4.1. Estimate of Confusion 

Previous surveys with Spitzer are not deep enough to 
measure the photometric confusion noise at 70/im, and the 
estimated confusion level is based on the observed bright 
source counts and models of galaxy evolution (Dole et al. 
2004b). With the ultra-deep data of GOODS-N, we can di- 
rectly measure the confusion level at 70/im. We define the 
instrument noise (including photon noise, detector noise, 
and noise associated with the data processing) as aj. The 
total noise (ctt) represents the noise after the extraction 
of sources above a limiting flux density (Sum), and the 
confusion noise (ctc) represents fluctuations due to sources 
with flux densities below Sum- As defined here, cjc is the 
"photometric" confusion noise, following the terminology 
of Dole, Lagache, & Puget (2003). In the direction of 
GOODS-N, the contribution of Galactic cirrus to the con- 
fusion noise is negligible (^ 0.01-0.02 mJy), based on the 
relationship given by Dole et al. (2003) and the calcula- 
tions of Jeong et al. (2005). 

The instrument noise was estimated empirically by sub- 
tracting pairs of data with the same integration time 
and covering the exact same region on the sky to re- 

5S70=S,.(71.4/im) and S24=S^(23.7/im) throughout this paper. 



move sources and any remaining residuals from the sky 
after filtering. The measured instrument noise integrates 
down nearly with t~°-^ (Fig. 1). For these data, we 
find fjf a i-i(l + l3t°-^), where /? = 0.04 for integra- 
tion time t in units of ks. The fimctional form of this 
relationship is based on empirical results from several dif- 
ferent data sets, and the /3 parameter depends on the 
background level and the quality of the data reduction. 
We use the above function to extrapolate the instrument 
noise from half the data to the full data set and derive 
(71 = 0.0399 ± 0.0036 MJysr-i. 

Since the total noise image (after source extraction) and 
instrumental noise image have nearly Gaussian distribu- 
tions, the confusion noise can also be approximated by 
a Gaussian and is given by Cc = ((t|, — aj)'^-^. We it- 
erate between source extraction at different limiting flux 
densities and confusion noise measurements until we con- 
verge to a solution with q = Sum/ac = 5. For the 
q = 5 solution, we derive or = 0.0485 ± 0.0034 MJy sr-^ 
and a,. = 0.0276 ± 0.0079 MJysr"^, for a limiting source 
flux density of S70= 1.5 mJy.^ Including the additional 
systematic uncertainties of the absolute calibration scale 
(10%) and the conversion between point source noise and 
surface brightness noise (10%, §3), we derive a point source 
confusion noise of ctc = 0.30 ± 0.15 mJy {q = 5). 

In comparison, the predictions of Dole et al. (2003, 
2004b) suggest a 5 = 5 photometric confusion level of 
(7c ~ 0.28 mJy, depending on the exact shape of the dif- 
ferential source counts. The measured confusion level of 
0.3 mJy agrees well with the predicted photometric con- 
fusion level. However, the source density criterion (SDC) 
confusion limit of 3.2 mJy adopted by Dole et al. is more 
than a factor of two higher than the limiting flux density 
derived here. The Dole et al. SDC limit corresponds to 
g ~ 7 and a high completeness level of > 90%. Sources can 
be extracted at lower completeness levels, and counts can 
be derived reliably well below this SDC limit by making 
use of the 24/im data to help extract the faintest sources. 

4.2. Source Counts 

The source counts were derived separately for the cen- 
tral 10' X 10' ultra-deep field, for the wide 0.614 deg^ GTO 
deep field, and for the 12'.9 x 12'.9 intermediate field which 
includes the ultra-deep field and the surrounding regions 
of intermediate depth between the GTO and ultra-deep 
surveys (Table 1). Analysis was done on the "intermedi- 
ate" field for better statistics at S70> 4 mJy. The central 
ultra-deep area has a deficiency (~ 30-50%) in number of 
sources with fiux densities between 5-10 mJy in compari- 
son to typical areas surrounding GOODS-N, presumably 
due to cosmic variance. The intermediate and GTO fields 
allow the measurement of source counts for flux densities 
(4-12 mJy) not weU sampled by the ultra-deep GOODS-N 
and the shallow xFLS surveys. The combination of the 
GOODS-N and the xFLS data sets yields source counts 
over the flux density range from 1.2 mJy to 455 mJy (Fig. 
2). In comparison, the counts presented here are 12.5 times 
deeper than the counts previously published by the MIPS 
Team (Dole et al. 2004a) and about 7 times deeper than 
the counts based on the xFLS verification field (Prayer et 
al. 2006). 
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Sources were extracted for signal-to-noise (S/N) ratios 
of S/N>2.5 for the ultra-deep field and S/N> 4 for the in- 
termediate and GTO-deep fields. The uncertainty image 
(§3) was used to represent the point-source noise as a func- 
tion of position. The number of spurious sources for each 
flux bin was estimated by performing source extraction on 
the negative image. Reliability degrades significantly be- 
low S/N< 4. To compensate for spurious sources at low 
S/N, we required the presence of 24/im counterparts in 
the ultra-deep field. Since the GOODS 24^m data (M. 
Dickinson et al., in preparation) are about 100 times more 
sensitive in flux density, very few, if any, 70/im-only detec- 
tions are expected. Within the ultra-deep field, the chance 
coincidence within 8" (~ 2 times the positional rms uncer- 
tainty of the faintest 70/im sources) of a 24/im counterpart 
with S24> 80/zJy is 20%. By requiring 24^m counterparts, 
we removed 80% of the spurious sources, and the observed 
counts were then corrected assuming a 20% chance match 
for the remaining spurious sources. 

After correcting for reliability, the observed source 
counts are corrected for completeness and flux biasing (i.e., 
the scattering of faint sources into brighter flux density 
bins) using the Monte Carlo approach described by Chary 
et al. (2004). Simulated point sources spanning the full 
range of flux densities were injected at random positions 
into each image (one at a time with 3x10'' repetitions per 
field). By using the simulations and the same source ex- 
traction methods, we populate a Pij matrix representing 
the probability that a galaxy with input flux density i is de- 
tected and measured with output flux density j. By sum- 
ming over the elements of the Pij matrix, the true counts 
are derived from the observed counts. This Monte Carlo 
approach allows the derivation of the counts at low com- 
pleteness levels; the completeness level for the faintest flux 
bin (1.2-1.6mJy) is only 30% (Table 1). The uncertainties 
on the counts include the uncertainties of the reliability 
corrections, the Poissonian errors propagated through the 
Pij matrix, and the 10% absolute flux calibration uncer- 
tainty combined in quadrature. 

The Euclidean-normalized differential source counts 
turn over around 8-10 mjy (Fig. 2). The observed counts 
are consistent within 30% of the Lagache et al. (2004) 
model for bright flux densities. We find a slightly larger 
number of faint galaxies (S70<3mJy) than predicted by 
Lagache et al. (2004), which may be due to cosmic vari- 
ance in the direction of GOODS-N. Ultra-deep 70/im ob- 
servations over larger areas and along different lines of 
sight are needed to constrain the models more accurately 
at faint flux densities. 

At low flux densities, a weighted least-squares fit to 
the differential source counts yields dN/dS oc with 
a = 1.6 ± 0.6. This is consistent with the power law of 
a = 1.6 ± 0.1 derived for the faint 24/im sources (Chary 
et al. 2004). The observed fiuctuations also constrain the 
faint source counts. We estimate the extrapolated source 
counts down to the confusion limit via simulations. Dif- 
ferent populations of sources covering a wide range of a 
values and normalizations at 2 mJy were randomly injected 
into the instrumental noise image and then extracted us- 
ing the same techniques carried out for the confusion mea- 
surement. The best-fit solution for simulations consistent 
with the observed constraints is a = 1.63 ± 0.34. At the 



confusion level of 0.3 mJy, we derive an extrapolated value 
for the Euclidean-normalized differential source counts of 
{dN/dS)S^-^ = 290 ± 200galsr-i Jy^-S. 

4.3. Extragalactic Background Light 

The expected extragalactic background light (EBL) due 
to infrared galaxies at 70/im is fairly uncertain due to 
the difficulty of deriving accurate zodiacal light correc- 
tions for the Diffuse Infrared Background Experiment 
(DIRBE) measurements (e.g., Wright 2004). By inter- 
polating between the 24/im EBL value (Papovich et al. 
2004), the 60/im value (Miville-Dechenes, Lagache, & 
Puget 2002), and the DIRBE 100/im and 140/im mea- 
surements (Wright 2004), we estimate a predicted EBL 
level of fliy = 10 ± 5nW m~^sr~^ at 71.4/im. In com- 
parison. Dole et al. (2006) estimates a total EBL of 

7.1 ± 1.0 nW m~^sr~^ from 70/im stacking analysis and 
the extrapolation of the 24/im counts, and the Lagache et 
al. (2004) model predicts a value of 6.4 nW m^^sr^^. 

By summing over the observed source counts for S70> 

1.2 mJy (including a small contribution from sources 
brighter than 455 mJy based on the Lagache et al. 2004 
model), we derive a contribution of 4.3 ± 0.7nW m~^sr~^ 
to the EBL. This is about 60% of the total EBL. By ex- 
trapolating {a = 1.63 ± 0.34) the 70/im source counts 
down to the confusion level, we derive a contribution of 
5.5 ± 1.1 nW m~^sr~^ for S70> 0.3mJy, and the extrapo- 
lation to zero flux density yields an estimated total EBL 
of 7.4 ± 1.9 nW m~^sr~^. The uncertainties on the EBL 
measurements include an additional 10% systematic un- 
certainty to the error budget, accounting for the uncer- 
tainties associated with the absolute calibration and color 
corrections. 

5. CONCLUDING REMARKS 

Based on ultra-deep 70/im observations, we derive 
source counts down to a flux density of 1.2 mJy, directly 
resolving about 60% of the EBL. The total fraction of the 
EBL estimated for sources down to the confusion level 
((Tc — 0.3 mJy, q = 5) is about 75%. A power-law ex- 
trapolation to zero flux density implies a total EBL of 
7.4 ± 1.9 nW m~^sr~-' at 71.4/im. This is consistent with 
the value predicted based on EBL measurements at other 
wavelengths, the value predicted from the Lagache et al. 
(2004) model, and the value derived from the extrapola- 
tion of the 24/im counts and 70/im stacking analysis (Dole 
et al. 2006). However, the uncertainties on the results 
leave open the possibility of a significant population of 
sources at low 70/im flux densities that are not accounted 
for in the models, such as highly obscured z 1 AGNs as 
proposed to account for the hard X-ray background (e.g., 
Worsley et al. 2005). Studies of the counterparts of the 
faint 70/im population are ongoing and will help to con- 
strain the infrared luminosities and the relative fraction 
of AGN versus starburst-dominated galaxies in the high- 
redshift Spitzer-selected surveys. 

We thank our colleagues associated with the Spitzer mis- 
sion who have made these observations possible. This work 
is based on observations made with the Spitzer Space Tele- 
scope, which is operated by the Jet Propulsion Laboratory, 
California Institute of Technology under NASA contract 
1407. 
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Fig. 1. — Noise for 4" pixels as a function of integration time (It). The corresponding point-source noise in mjy is shown at the right. 
Instrumental noise measurements are shown as open circles and axe represented by the solid line. For comparison, the dotted line shows a 
f-0.5 fmiction. The derived confusion level (uc) is shown by the dashed-dotted line, and the total noise after the extraction of sources with 
570 > 5(7c is shown by the dashed line. The total noise and instrument noise for the ultra-deep field are shown by the solid square and solid 
circle, respectively. 
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Fig. 2. — The 70fim differential source counts from the ultra-deep GOODS-N field are shown as filled circles, the intermediate field counts 
are shown as open triangles, while the results from the GTO deep field are shown as open circles. The xFLS number counts (Fraycr et 
al. 2006), corrected for the updated calibration, are shown as squares. The vertical bars represent the total errors, and the horizontal line 
segments show the sizes of each bin (Table 1). The model of Lagache et al.(2004) is shown by the solid line. The dashcd-line represents 
the best fit extrapolation of the counts {dN/dS oc S~'^) down to the confusion level shown by the solid triangle, and the allowable range of 
extrapolated counts constrained by fluctuation analysis is given by the dotted-lines. 
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Table 1 

GOODS-N 70fim Source Counts 



S(low) 
[mJy] 


S(high) 
[mJy] 


S(avg) 
[mJy] 


Observed 
Number 


Reliability 


Completeness 


Corrected 
Number 


Log[{dN/dS) 5-2-5] 
[gal sr~^ Jy^'5] 


1.2 


1.6 


1.38 


21 


0.75 


0.30 


52.4±14.3 


3.043±0.136 


1.6 


2.2 


1.87 


28 


0.85 


0.49 


48.7±11.0 


3.163±0.118 


2.2 


3.0 


2.49 


31 


0.97 


0.75 


40.2± 7.8 


3.267±0.107 


3.0 


4.5 


3.60 


21 


0.97 


0.80 


25.6± 6.1 


3.197±0.122 


4.5 


8.0 


6.29 


14 


1.0 


1.05 


13. 3± 3.6 


3.150±0.133 


8.0 


15.0 


11.68 


11 


1.0 


1.04 


10.6± 3.2 


3.423±0.147 


2.2 


3.0 


2.55 


28 


0.75 


0.43 


48.7±13.4 


3.152±0.136 


3.0 


4.5 


3.73 


26 


0.69 


0.54 


33.5±10.8 


3.129±0.154 


4.5 


6.5 


5.73 


20 


1.0 


0.86 


23. 2± 5.2 


3.311±0.117 


6.5 


11.0 


8.81 


19 


1.0 


1.10 


17.3± 4.0 


3.300±0.119 


11.0 


16.0 


12.62 


10 


1.0 


0.87 


11.5± 3.7 


3.466±0.153 


6.0 


7.5 


6.68 


101 


0.86 


0.54 


160.1±21.9 


3.318±0.088 


7.5 


9.0 


8.19 


96 


0.95 


0.69 


132.5±16.8 


3.458±0.085 


9.0 


11.0 


9.89 


73 


0.99 


0.71 


101.4±13.2 


3.421±0.086 


11.0 


14.0 


12.07 


62 


1.0 


0.92 


67.3± 8.7 


3.283±0.086 


14.0 


18.0 


15.72 


45 


1.0 


0.95 


47.5± 7.4 


3.294±0.094 


18.0 


30.0 


23.48 


43 


1.0 


1.09 


39.6± 6.1 


3.173±0.094 



Note. — The first 6 rows are for the ultra-deep field (0.0277 deg^), the next 5 rows are for the intermediate field 
(0.0462deg2), and the remaining rows are for GTO-deep field (0.614deg2). The corrected counts are equal to the observed 
counts multiplied by the reliability value and divided by the completeness correction. The uncertainties for the Euclidean- 
normalized differential counts [{dN/dS) S'^ '"'] include the Poissonian noise, the uncertainties associated with the reliability, 
completeness, and flux biasing corrections, and the absolute calibration uncertainty of 10%. 



